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A random matrix theory approach is applied in order to 
analyze the localization properties of local spectral density for 
a generic system of coupled quantum states with strong static 
imperfection in the unperturbed energy levels. The system is 
excited by an external periodic field, the temporal profile of 
which is close to monochromatic one. The shape of local 
spectral density is shown to be well described by the contour 
obtained from a relevant model of periodically driven two- 
states system with irreversible losses to an external thermal 
bath. The shape width and the inverse participation ratio 
are determined as functions both of the Rabi frequency and 
of parameters specifying the localization effect for our system 
in the absence of external field. 



I. INTRODUCTION 

The statistical properties of many-body quantum ob- 
jects attract considerable attention in a broad range of 
modern physics including condensed matter physics and 
quantum optics. Of a special interest are the proper- 
ties specifying the temporal evolution for mesoscopic sys- 
tems of quantum states coupled by interaction. Extensive 
investigations of many-body interacting systems (such 
as nuclei, many-electron atoms, quantum dots, quan- 
tum spin glasses and quantum computers) have shown 
that the state-state interaction can be responsible for 
the quantum localization effect and the quantum chaos 
in isolated complex systems [1-9] . As in Random Matrix 
Theory (RMT) [10,11], the level spacing statistics is well 
described by the Wigner-Dyson distribution and the in- 
dividual basis states are spread over the large number of 
eigenstates. In a sense the interaction leads to dynamic 
thermalization without coupling to an external thermal 
bath. On the other hand, to our knowledge, there are no 
direct theoretical studies for mesoscopic systems driven 
by an external field. 

In the work we apply a RMT approach in order to 
investigate the localization effect for a complex system 
of coupled quantum states in the presence of an exter- 
nal field. In the approach a single matrix from a given 
statistical ensemble (characterized by a few matrix pa- 
rameters) represents the Hamiltonian H (t) for the typ- 
ical (generic) dynamical system of given class. Our sys- 
tem is represented by a matrix, which includes a lead- 
ing diagonal with disordered random values and random 
off-diagonal elements inside a band with the halfwidth 



b. This band random matrix with disordered diagonal 
(BRMDD) is a generic Hamiltonian model for the sys- 
tems with strong static imperfection in the unperturbed 
energy levels. The imperfection plays an essential role 
for the quibit composition and can destroy the operabil- 
ity of quantum computer [9,12,13]. BRMDD-based mod- 
els were applied to study the electron transport problem 
[14-16] and the problem of interacting particles in a ran- 
dom potential [17-19]. The results obtained with the help 
of these models are of an obvious interest also in analyz- 
ing such few-freedoms physical objects as the vibrational 
quasicontinuum of polyatomic molecule [20-22]. 

The statistical properties of conservative BRMDD- 
based systems were studied in some detail [23-29]. The 
shape, the localization length and the inverse partici- 
pation ratio (IPR) for eigenfunction have been investi- 
gated with the help of numerical simulations [23,27,28] 
and the supersymmetry approach [28,29]. These investi- 
gations have exhibited the Lorentzian shape for the local 
spectral density (LSD) of states under the circumstances 
when a non-perturbative localization regime is realized. 
The state-state interaction strength, at which the eigen- 
states are extended over the whole matrix size N and the 
eigenenergy level spacing statistics has the Wigner-Dyson 
form, have been revealed for the BRMDD with essentially 
large band (when 25 -I- 1 > y/N) [26-29]. The associa- 
tion between ergodic properties of LSD and the matrix 
parameters have been clarified recently for the BRMDD 
with arbitrary small band [30]. The lack of ergodicity 
for LSD is shown in [30] to be identified with an expo- 
nential increase in IPR with the strength of state-state 
interaction. 

We study the LSD of states for a generic BRMDD- 
based quantum system to be excited by an external quasi- 
monochromatic field. The LSD was introduced in 1955 
by Wigner [31,32] and successfully employed in RMT to 
describe statistically the localization effects for complex 
quantum systems [10,11] (including the systems repre- 
sented by band random matrices with reordered leading 
diagonal [33,34]). In our study this quantity specifies 
spreading of the energy concentrated initially in an in- 
dividual basis state \g), between the quasienergies and 
gives the quasienergy spectrum of \C {lu)\ where C (uj) 
is the Fourier transform for correlation 

Cit)^{g\expi~iU{t)t/h)\g). (1) 

As for the conservative systems [10,11], the shape of LSD 
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may be characterized by the width F of a quasienergy 
scale, on which the individual state \g) is localized. The 
number of quasienergies populating this scale is given 
by the ratio F/A^^ = L designated here as a localiza- 
tion length of LSD (A^j is the mean quasienergy spac- 
ing). Hence the quantity L specifies the greatest possible 
number of quasienergies, where the basis state |g) ca: 
be effectively admixed. The ergodic properties for th 
system can be identified with the structure of LSD. Th 
non-crgodic LSD is a strongly fluctuating spiked func 
tion and the IPR ^, which gives the actual number c 
quasienergies involving the state \g), is low in compari 
son with L. In the crgodicity case the LSD is monotoni 
and the number ^ approaches the value of L. 

Here we study the shape and the ergodic features c 
LSD and estimate the shape width F and the IPR 
as functions of field-system interaction strength. Th 
functions are analyzed in relation to relevant propertie 
(the mean quasienergy spacing, the width, the localize 
tion length and the IPR) specifying the localization effec 
for our BRMDD-based system in the absence of externg 
field. 



II. MODEL DESRIPTION 

We consider a generic system of quantum states, th 
unperturbed energy levels of those are depicted in the 
inset of figure 1. The system consists of a single state \g) 
andN = 2K+1 states |A;) {k = -K, ...,K). The total 
Hamiltonian of the system is considered as the sum 



H {t) = H(°) + H(") (t) . 



(2) 



Here the time-independent operator H^^^ describes a 

state-state interaction in the absense of field. The time- 
dependent part H'^^) (t) specifies a field-system interac- 
tion. 

In the basis of unperturbed states the Hamiltonian 
jj(o) represented in terms of a real symmetric matrix 
with statistically independent random elements 



I H(") \j) 



+ ^ 



V|j) 



(3) 



where the off-diagonal elements («|V|j) 
specify the state-state interaction = g, —K, . . . , K). 
In the model we take into account an interaction between 
the states \k). The corresponding elements Vfc'fc are dis- 
tributed uniformly in the interval [—V, V] with (Vfc'fc) = 
and (V"fc,fc) = = v'^ if \k' — k\ < b or are zero other- 

wise. The coupling of the state \g) with \k) is ignored and 

Vgk = Vkg = 0. The diagonal elements E^^^ correspond- 
ing to the energy levels of states \k) are uniformly dis- 
tributed according to the Poisson statistics with the mean 
spacing A between adjacent levels: —KA < E^'' < KA. 
One of the energy levels (for definiteness, we take the 



level Eq^'^ for the state |0)) is located in the midpoint 

= 0. The state 



,(0) 



of interval [-KA,-KA]: (^E^ 
\g) is the lowest one. The energy levels are chosen so that 



-^o'^ - E^"' > 2KA. Note that in the basis of states \k) 
the operator H^^^ is represented by a BRMDD. 
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FIG. 1. A temporal shape of external periodical field (solid 
lines, M = 16). The fat dashed line gives the envelope 

sin{u! ft). The inset shows a diagram of unperturbed energy 
levels for represented BRMDD-based quantum system. 

An external periodic field with the frequency cuj = 
Eq'^ — Eg^^ stimulates the transitions between the non- 
perturbed basis states l^) and |0). The field-induced cou- 
pling between \g) and the states |fc) with fc 7^ is ac- 
cepted to be negligible. The temporal profile of the field 
is shown in figure 1 to be a piecewise function with the 
envelope sin^LUft). The total number of pieces covering 
the field period interval T/ = 27r/a;/ is equal to M. At 
large M the field shape is close to monochromatic one. 
On each of the time pieces [tm, im+i] (the integer m gives 
the number of piece interval) the field is time-invariant 
and the field-system interaction is represented by 



H(") (tm) = Osin(27rm/M)D. 



(4) 



Here the Rabi frequency O specifies the strength of field- 
system interaction. The operator D specifies the field- 
induced coupling between the states: 



(^1 D \j) = S^Sgj + SigSoj. 



(5) 



At such an approximation the evolution of system for one 
field period interval is described by a unitary operator 



M 



U(T/) = n exp {-i [H (tm)] Tf/Mh) . 



(6) 
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From diagonalization of the matrix \J{Tf) we obtain the 
eigenf unctions \aj) and the quasienergies wj. The quan- 
tity Wjg = I {cxj\g) P gives a probabiUty to find the prob- 
ing basis state \g) in the eigenstatc \aj). 

The described model can be considered as a driven two- 
state system (the states \g) and |0)), the upper state |0) 
of which is coupled by interaction with a number of oth- 
ers quantum states. The interaction-induced probability 
distribution of the state |0) over the eigenstates of H^^^ 
corresponds to a quasicncrgctic distribution of |0) in the 
absence of field (fi = 0). Since the operator H^"^ is given 
in terms of a BRMDD the probability distribution of |0) 
over quasienergy scale is described by the Breight-Wigner 
shape (that is, by the Lorentzian contour) [23,27,28,30] 



Po (w) = 



1 



277 w2 + rl/4 



(7) 



of the width Fq. In some sense the state |0) is Lorentzian- 
broadened. Therefore we associate our BRMDD-based 
system with a two-state system, the upper state of which 
is continuously Lorentzian-broadened with the width Fq 
because of irreversible losses to an external thermal bath. 
Here the lower state has no losses and hence is not 
broadened. This two-state system with irreversible losses 
(TSSIL) is stimulated by an external monochromatic 
field. Such a well-known model (sec, for instance, [35,36]) 
corresponds to our BRMDD-based system in the limit of 
infinitely large number N of coupled states \k). 

Depending on the ratio between the width Fq and the 
Rabi frequency f2 the TSSIL model gives two shapes of 
\C [35,36]. At 20/ro < 1 the correlation C{t) is an 
decaying aperiodical function of time. Then the shape of 
|C (w)l^ is described by a contour 



\C{u;)\ 
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{uj^ + al/4) (uj^ + all 4) 



(8) 



with coefficients a\<^ = F§/2 - ± V^^f^y^^^W). 
At a significantly weak field (2ri/Fo <C 1) this shape is 
close to the Lorentzian contour, the width F of that is a 
quadratic function of the Rabi frequency f2: 



2f2Vro 



(9) 



That implies a single exponential decay with time for 
C{t). At a stronger field (2i^/Fo > 1) the shape of 
|C (w)l^ is given by 



|C(a;)|^ = ^ 



(4a;2 



dl 



(10) 



where df = fl"^ -Tq/A and (i2 = ro/2. Then the correla- 
tion C{t) shows a decaying oscillation. In the limit of a 
strong field (20/ro ^ 1) the shape width F is invariant 
to the Rabi frequency and approximates the magnitude 

ro/2. 



III. NUMERICAL SIMULATION 



The properties of LSD are analized for a wide range 
of parameters N, v, A and Q {3 < N < 1000, 10"^ < 
v/A < 10^ and lO'^ < n/A < lO^). The total num- 
ber M of time pieces covering the period interval T/ is 
32. Tli(^ shape and the IPR of LSD are determined from 
averaging over disorder (that is, over many random ma- 
trices). The number of disorder realization is more than 
200. 



A. Shape of LSD 

We restrict the study by the case of moderate strengths 
of state-state interaction when the width F is essen- 
tially small compared to the whole quasienergy scale 
(F <C NA/U). It implies that the perturbation- induced 
variations in level density are negligible and the quasinen- 
ergies w„ are homogeneously distributed in the energy 
band [—KA^^jKA^^] with the mean quasienergy spacing 
= A/Ti. Then the shape of LSD can be defined as 



p{uj) = 



1 {i:,w,,5{^-. 



x;,'5(w-wi) 



(11) 



where (...) means the averaging over disorder. We asso- 
ciate the shape p (w) with contours represented in equa- 
tions (8) and (10) and determine the coefficients ai,2 and 
di 2 by mean square fitting of equation (8) or (10) to an 
averaged LSD. The shape width F is estimated from the 
obtained coefficients as a twice halfwidth for the right 
slope of p {lij) in the positive part of quasienergy scale. 

The numerical simulation for different magnitudes of 
A^, Fq, ^0 and VL shows that the properties of p{u)) 
are determined by some factors. First of all, that is 
the mean quasienergy spacing A^^ . Another important 
factor is the width Fq of shape po (w) for a probability 
distribution of |0) over quasienergies in the absence of 
field. We should distinguish two localization regimes to 
be realized for the distribution. A perturbative local- 
ization regime takes place when the state-state interac- 
tion is weak and the state |0) is concentrated mainly in 
an individual quasienergy. For this regime the distribu- 
tion width Fq is small in comparison with A^^. On the 
other hand, the strong state-state coupling results in a 
broad spreading of |0) over the quasienergy scale. Such 
a regime (designated here as a non-perturbative one) is 
characterized by large magnitudes of Fq: Fq > A;^. Cer- 
tainly, the third significant factor charachterizing the fea- 
tures of p {lo) is the strength Vl of field-system interaction. 
Our study demonstrates that the ergodic properties to be 
specified in terms of IPR ^0 for the quasienergetic distri- 
bution of |0) have no effect on the shape p {uj). Therefore, 
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wc consider the shape width F as a function offl, Fq and 
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FIG. 2. Shape p{uj) of LSD at 2n/r(, = 0.5 (Q), 2(A), 
4(V). The hnes show shapes obtained from fitting of contours 
(8) (dashed line) and (10) (solid lines) to p{uj). 

Similar to the TSSIL model, at a weak field {2^1 /Tq < 
1) the shape is shown in Figure 2 to be well de- 

scribed by contour (8). The obtained values of P sat- 
isfy the requirement ri^/Fo < F < 20. The type of 
f2-depcndcncc for the width F is specified by relation be- 
tween Ai^, Fo and ft. Our study demonstrates that the 
low values of Fq or Q give the linear low for the depen- 
dence. At Fo < the width F is seen from the inset of 
Figure 3 to be a linear Fo-independent function of O: 



F w 2f2. 



(12) 



This dependence can be considered as a manifestation 
of the perturbative regime for distribution po {^)- At 
Fq > (when a non-perturbative localization regime 
for po {^) is realized) the increase of O results in the 
transition of 0-dependence for F from the linear law to 
the quadratic one. We associate such a behavior of the 
width F with the magnitude of localization length L for 
LSD. Figure 3 shows that a linear function of O 



T^An^/A^ 



(13) 



takes place when L <C 1. The coeflacient A — 0.7 is ob- 
tained from mean square fitting of (13) to calculated val- 
ues of F. On the other hand, at L » 1 the f2-dependence 
for the quantity F is close to the quadratic low (9) pre- 
dicted by the TSSIL model. As for distribution po (w), 
the transition from the linear to quadratic f2-dependence 
of F can be associated with crossover of the shape p (w) 
from the perturbative to non-perturbative localization 
regime. 




n/(r A )' 

FIG. 3. The localization length L as a function of 
n/VA^Fo at a weak field (20/Fo < 1) for Fo/A<^ = 0.9 
(squares), 2.5 (V), 20 (stars) and 225 (full circles). The 
dashed line gives approximation (13) obtained for the points 
with Fo/Au, > 1. The solid line represents law (9) predicted 
by the TSSIL model. The inset shows the ratio F/Fo as a 
function of 2n/Fo at 2n/Fo < 1 for To/A^ = 0.06 {+), 0.2 
(O), 2.5 (A), 20. (stars) and 226 (full circles). The dashed 
line corresponds to dependence (12). The solid line represents 
law (9) predicted by the TSSIL model. 
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FIG. 4. The ratio F/Fo as a function of 2n/Fo at a strong 
field (2n/F(, > 1) for Fo/A^ > 1 (full circles) and Fo/A^ < 1 
(O)- The solid line shows a relevant dependence predicted by 
the TSSIL model. 
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At a strong field (2ri/ro > 1) the shape p {oj) is shown 
in Figure 2 to be well fitted by contour (10) obtained 
from the TSSIL model. As for the case of weak field, 
the behavior of ^-dependence for the shape width T is 
determined by the localization regime realized for shape 
po (t^)- Figure 4 demonstrates that at Fq > A^j the O- 
dependcncc of T is in good agreement with a relevant 
curve obtained from the TSSIL model. At high field am- 
plitudes (when 20/Fo 1) the width F approximates 
the value ro/2. The minor discrepancy between T and 
the curve predicted by the TSSIL model can be explained 
by a divergence of the temporal profile of field from 
the monochromatic shape. At the perturbative regime 
of po (w) the f2-dependence of F is proportional to the 
TSSIL curve. Our analysis shows that for any magni- 
tude of Vt the calculated values of F are higher then ones 
from the TSSIL model by a factor of 1.4. In the limit of 
high field amplitudes the width F is close to the magni- 
tude Fo/1.4. 



B. Inverse psirticipation ratio of LSD 

Wo investigate the IPR ^ under the circumstances of 
localized regime for the probability distribution of basis 
state |0) over the eigenstatcs of H^"'. It implies that the 
state |0) is spread over a sufficiently large number of the 
eigenstates and ^ ^ 1- Here the quantity is the 
IPR for the probability distribution of |0) and gives the 
actual number of eigenstates involving the state |0) [30]. 
In the study the main attention is paid to the case of 

C » 1. We accept the IPR as ^ = ((Ej and 
obtain this quantity from numerical simulation. Here 
(. . .) means the averaging over disorder. The IPR ^ is 
considered as a function of f7, and Fq. Results of the 
study are represented in Figure 5. 

Our study demonstrates that at fixed O and Fq the IPR 
^ varies proportionally with the parameter ^o- This fact 
testifies that the change in ^ with the Rabi frequency is 
due to transformations in the shape of LSD. For instance, 
in the limit of strong field {2Vl/Tq > 1) the width P of the 
shape p {lo) is fl- independent and F « Fo/2. As a result, 
the quantity ^ is also invariant to Q. and approximates the 
magnitude of (see Figure 5). In other words, the IPR 
of LSD is equal to the IPR for the distribution of |0) over 
the eigenstates of H^^^ (here wo should take into account 
the presence of two peaks of LSD in the negative and 
positive parts of quasienergy spectrum). As indicated 
earlier, at a weaker field (2f7/ro 1^ 1) the variation in 
give rise to essential shape transformations both in the 
shape width and in the contour form. The analysis of 
calculated data shows that the behavior of IPR for this 
parameter region can be empirically approximated by the 
law: 



where the coefficients B = 1.54 and /? = 1.18 are ob- 
tained from mean square fitting of (14) to calculated val- 
ues of ^. 
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FIG. 5. Dependence of the IPR ^ on the ratio 2Q/Tq at 
Po/A^ = 40 (Co = 30 (+), 40(stars)) and Vq/A^ = 100 
(Co = 39 (0)i 100 (full circles)). The dashed line shows ap- 
proximation (14). 

Notice that the observed Jl-dependence of IPR is close 
to the linear law and differs significantly from the ex- 
pected quadratic dependence to be realized in the case of 
the Lorentzian shape for LSD. Probably, this difference 
is attributed to the fact that we are not far enough in the 
asymptotic regime of weak field implying 20, /Fq ^ 1. In 
our simulation we have 2$7/ro ^ 0.1 and the shape of 
LSD appers to can not be well reduced to the Lorentzian 
contour. Unfortunately, the numerical simulation of lo- 
calized regime for LSD with 2il/ro < 0.1 requires too 
large matrix sizes and much computational efforts being 
beyond our numerical abilities. 



IV. CONCLUSION 

We have analyzed localization and ergodic properties 
of LSD for a periodically excited generic system of cou- 
pled quantum states with strong imperfection in the un- 
perturbed energies. These properties have been demon- 
strated to be determined essentially by a state-state in- 
teraction resulting in the localization effect for the system 
in the absence of external field. If the interaction is so 
strong that the conservative system exhibits this effect 
the shape of LSD can be obtained from a relevant model 
of driven two-state system with irreversible losses. In a 
sense the state-state interaction may be recognized as the 
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losses acting till the mcsoscopic effect becomes evident. 
The time scale and the rate of a correlation decay caused 
by such losses are specified in terms of the shape width 
and the IPR of LSD. 

Author thanks V. Churakov for useful discussions and 
remarks. 
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